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Abstract : In this paper we consider a viscoelastic three dimensional body 
(of Maxwell-Boltzmann type) controlled on (part of) the boundary. We assume 
that the material is isotropic and homogeneous. If the body is elastic (i.e. no 
dissipation due to past memory), controllability has been studied by several au¬ 
thors. We prove that the viscoelastic body inherits the controllability properties 
of the corresponding purely elastic system. The proof relays on cosine operator 
methods combined with moment theory. 


1 Introduction 


We consider a linear viscoelastic body occupying a region P C ]R^. We assume 
that the body is isotropic and homogeneous so that the dynamics of the body 
is described by the operator £: 


£u = ^Au -I- (A -I- ^)V (V • u) , u = u(a:) x G P . 

Here u G IR^ and \ > 0, fj, > 0 are the Lame coefficients. We assume that A 
and fj, are constant (i.e. that the body is homogeneous). 

Boldface denotes vectors. Note that also the space variable (i.e. x) is a 
vector, with dimx = dimu, but we don’t use boldface for the space variable 
(and not for the vector 0). 

If the body is elastic, the evolution in time of the displacement is described 
by the Navier equation 


u" = /lu -I- F , 


u(0) = uo 
u'(0) = ui 


( 1 ) 


(u = u(x,t), F = F(x,t) and x G H C IR^) while if the body is viscoelastic 
(of the Maxwell-Boltzmann type) the evolution in time of the displacement 
(denoted w = w(x,f) G H^) is described by 

w” = Cw + M{t- s)£w(s) ds -t- F , | (2) 
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(see [mp. 112-113]) In fact, the integral should extend from — oo but in the 
study of reachability for linear systems, the control being applied after a time 
to, say to = 0, it is not restrictive to assume w = 0 for t < to = 0- 

We apply a displacement on a part F of the boundary of fl (F = dil is 
not excluded), both to the elastic and the viscoelastic body, and we are go¬ 
ing to study whether is it possible to control the pair {{velocity of displace¬ 
ment) / {displacement)) to hit a prescribed target at a certain time T. Fe. we 
impose the following boundary condition for t > 0: 

u(a;, t) = f (cc, t) a: € F u(a:, t) = 0 x £ dfl \ F (3) 

•w{x,t) = i{x,t) X £T ■w{x,t) = 0 X £ dfl\T (4) 

(if F = did then disregard the condition on dil \ F). 

When the initial condition and the affine term F are zero and we want to 
stress dependence on the control f, we write w^. 

Finder the assumptions we shall state below, it turns out that {u'{t),u{t)) 
and (’w'(t), w(t)) are [iF“^(fl)] ^ x [L^(n)] ^-valued continuous functions so that 
their values at a fixed time T make sense in this space. 

If the initial conditions are zero and also F = 0, the reachable sets at time 
T are the sets 

i?£;(T) = {((uf)'(T),uf(r)), f ^ 

Rv{T) = {{{w^y {T),w^{T)), f£L^ (0,T;[L2(r)]^)} . 

Note that (both in the purely elastic and in the viscoelastic case) the reachable 
set increases with time. 

Controllability at time T is the property Rv{T) = [iF“^(f2)]^ x [L^(n)]^ 

for the viscoelastic system, and Re{T) = [iF“^(n)]^ x [L^(r2)]^ in the elastic 
case. 

Under the assumptions we describe below, controllability in the purely elastic 
case, i.e. for the system m-m, has been studied by several authors (see the 
references in Section fL^ . We are going to prove that the controllability property 
which holds in the purely elastic case is inherited by the viscoelastic system. 

1.1 Notations, assumptions and the main results of this 
paper 

In this section we state the assumptions, we describe preliminary results on the 
controllability of the elastic system and we state our main results, which will be 
proved in the next sections. 

Assumption A) We assume: 

• the kernel M{t) is of class iliQ(.([0, -|-cxd)). 

• the region il C ]R^ is bounded and dil is of class . 

• We assume that both the Lame constants A and /i are positive and constant 
(i.e. the body is homogeneous). 

• The subset F of dil will be called the active part of the boundary. The 
first assumption on F is that it is relatively open in dil. 
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Known facts for the elastic system o with boundary conditions ([3]), prved 
in the references cited in Sect. O 

• Let T > 0 and f = 0. Let F € (o,T; uq e 

ui e Then 

t^(u(t),u'(t))GC([ 0 ,T];[ 7 Ji(f])] 3 ) xC([ 0 ,r];[L 2 (f])] 3 ) . 

• Let T > 0. Let F G (O, T; [^^(ri)]^), / G (q, T; [^^(r)]^) and Uq G 
[L2(f7)]3, ui G [H-\n)]^. Then 

t^(u(t),u'(t))GC([0,T];[L2(r!)]3) xC([0,T];[iJ-\L!)]3) . 

• There exist open subsets T of dVl and times T such that the purely elastic 
system is controllable in time T, i.e. 

RE{T)=[H-\n)f . (5) 

Assumption B) The active part T of dTl and T > 0 are chosen in such a way 
that the purely elastic system CD with boundary control ([3]) is controllable in 
time T, i.e. equality ([5]) holds. 

Remark 1 (On the notations) From now on, for the sake of readability, we 
drop the exponent 3 and we write simply L^(n), iLg(n), instead of 

[L^n)]\ 

We recall the following integration by parts formula, which we shall repeat¬ 
edly use: 


(£u) ■ (p dx = / (Tu) • ^ dr — / (T</>) • u dr -b / u • {Cep) dx . (6) 


The boundary operator 7~, the boundary traction, is: 


T0 = (n • V) 0 -b (A -b (V • u) n 


where n is the exterior normal to dCl. 

The following property, known as (output) admissibility or direct inequality 
is proved for example in Ch. IV]: 

Lemma 2 Let u solve m and m with f = 0. For every T > 0 there exists a 
number M such that 

|Tu|i 2 (o_T;L 2 (r)) < M -b |ui|^ 2 (Q) + |F|Li(o,T;L 2 (n))) ■ (V 

Note that the inequality usually proved is 

l7iu|L2(o_T;L2(r)) < M -b |ui|^2(n) -b |F|ii(o_-r.i2(f2))) 

(71 is the normal derivative on 9n) from which ([71) follows because u = 0 on 
dfl implies Vuk = 71 for every component Uk of u. 

Our first, ancillary result is as follows. 
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Theorem 3 Let Assumption A) hold. Then: 

1. Let F S L^(0,T;L^(r2)) and f = 0. System ^ with initial conditions 
w(0 ) = wq G HQ{n), w'(0) = Wi € L^(r2) admits a unique solution 
w G C'([0,r];iJoi(rj)) nc^ ([0,r];L2(rj)) for every T > 0. 

2. Let F G L^(0,T] L^{n)). System (0) with initial conditions w(0) = wq G 

L^{il), w'(0) = wi G and boundary control f G L^(0, T; L^(r)) 

admits a unique solution w G C ([0, T]; L^(r2)) n ([0, T]; for 

every T > 0. 

3. Let f = 0, Wo G Hq{LI), wi G Lp‘{Tl) and let w solve and Then, 
for every T > 0 and every T C dVL there exists M such that 


|Tw|^ 2 (o_r.i 2 (r)) < M (|wo|^i(q) + |wi|i 2 (Q) + |F|ii(o,T;L 2 (n))) • ( 8 ) 


This inequality is the direct inequalities of Eq. (0). 

We noted already that the statement of Theorem [3] holds when M = 0. 

The control result which we intend to prove is: 

Theorem 4 Let Assumptions A) and B) hold so that Re(T) = H~^(n)xL^{n) 
and let Ti > T. Then, system m, ® is controllable in time Ti, i.e. Rv{Ti) = 
H-^{Vt) X L‘^{Vt). 

The proof is based on two main steps. Step 1) proves that if the purely 
elastic system (HD is controllable at time T then the reachable set RyiT) of the 
viscoelastic system is closed with finite codimension in H~^{Ll) x L‘^{Ll). Thanks 
to the fact that the reachable set is increasing, this holds at every Ti > T. Step 
2) proves that if the purely elastic system is controllable at time T and Ti > T 
then we have 


Rv{Ti) = H-^{Tt) X L^{Tt). 


The organization of the paper is as follows: the next subsection recalls previ¬ 
ous results while section [5] presents preliminaries on the cosine operator theory 
and the proof of Theorem [31 

The two steps of the proof of Theorem [4] are in Section [d] 

1.2 References to previous work and preliminaries 

Controllability of the Navier equation (HD has been studied in several papers, of 
which we cite laiiiiiis]. The proofs in these papers are based on the inverse 
inequality, obtained using multiplier methods, i.e. it is proved that if F and T 
are suitably chosen then there exists m > 0 such that 



Here it is assumed f = 0, F = 0 and of course wq G wi G L^(fl). As 

noted in [31 p. 228], this inequality implies 



( 10 ) 
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(for a different m > 0) and this is the “dual” formulation of the definition of 
controllability. 

We mention that the inverse inequality in the paper is proved also if the 
Lame coefficients are (slowly) space varying and that isotropy is not assumed 
in [19]. These results can be extended to the viscoelastic case, precisely as we 
do below in the isotropic homogeneous case. 

In this paper, systems m and @ are studied using a cosine operator ap¬ 
proach and also moment methods. Let us define the following operators in 

domA = , A4, = CcI), A = i{-Af/^. (11) 


The operator A is selfadjoint with compact resolvent. Hence, admits an 

orthonormal basis of eigenvectors of A. Let it be denoted {</>„} and let — 
be the eigenvalue of 4>^ (it is known that the eigenvalues are negative and we 
assume that they are ordered in such a way that A„ < A„+i). The operator A 
generates a Co-group of operators on IR so that we can define the operators 


(the operator R+{t) is called the cosine operator generated by A). Note that 
t I—!> i? 4 .(t)u and t R-{t)u belong to C (L^(H), C([0, T], L^(r2))) and also to 

The operators R+{t) and R-{t) have the following expansions in series of 


^ +00 


+ 00 


R+it) I ^OLn(f)n I = X! (anCOsAnt) 0„(a;) , 

\n^l / n^l 

/ +CXD \ / +00 

R-{t) i'^an4>n] = * X! 


( 12 ) 


\n—l 


\n—l 


The solutions of problem (jH)-® are given by (see [2]) 


u(f) = i?+(<)uo-l-^ ^i?_(<)ui-|-.4 ^ f R-{t—s)F{s) ds—A f R-{t—s)Df{s)ds. 

Jo Jo 


(13) 


The operator D in this formula is the operator f u = Di where u solves 


f u(a;, t) = f{x,t) a: e r , 
u(a;, t) = 0 X G on \ T. 


It is known (see jU Th. 3.6]) that D G C (L^(9H), L'^iyi)) and, as noted in 
p. 796], it takes values in hence it is a compact operator. 

We use repeatedly the fact that when w G dom A we have 


7~w = —D* Aw (14) 

(a fact that can be proved using the integration by parts m as in El p. 181]). 

The cosine operator approach to controllability of systems with persistent 
memory, when u = rtGll, w = r(;Gll and C = A, was first used in [12] where 
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the existence of the control time was proved, but the control time itself was not 
identified. The control time was identified in dimension c? = 1 in subsequent 
papers and using moment methods, see for example m [H m [H and for the 
scalar valued wave equation in 17 C ]R^ in m- In this paper we combine the 
cosine operator method and the moment method in order to get a proof of 
controllability for Eq. This is based on the use of the following known 

estimates for the eigenvalues of the operator A (see BM)- there exist m > 0 
and M such that 

m < M ■ (15) 

2 Cosine operators and the proof of Theorem [3] 

Now we prove Theorem |3] in the case M ^ 0. The proof of the items [T] and B 
is similar to that in [m Chapt. 2] and it is only sketched. The first step is 
a definition of the solutions of problem (l2|)-(|4l). We first apply formally the 
MacCamy trick. Let R{t) be the resolvent kernel of M{t), given by 

R{t) + [ M{t — s)R{s) ds = M{t) . (16) 

Jo 

We “solve” the Volterra integral equation (2)) in the “unknown” Cw. We get 

Cw(t) = w"(t) — F(t) — f R{t — s)[’w"(s) — F(s)] ds . 

Jo 

We integrate by parts and we get 


w" = £w + aw'(t) + 6w(t) + f K{t — s)w{s) ds + Fi(t), 

t 

Fi(t) = F(t) + f i?(t — s)F(s) ds — i?(t)wi — i?'(t)wo , 

Jo 


a = R{0), b=R'{0), K{t)=R"{t) 


and the initial and boundary conditions (|3]). i.e. MacCamy trick removes the 
differential operator from the memory term. 


Remark 5 The term aw'(t) can be removed from the right hand side if we 
perform the transformation v(t) = e“^“/^^*w(t). The effect on the initial con¬ 
ditions is that v(0 ) = vq = wq while v'(0) = Vi = Wi — (a/2)wo and the 
control is replaced by (t). Of course this has no influence on control¬ 

lability and so with obvious modifications of the definitions of f(t), Fi(t), b 
and K{t) we shall work with (1171) but with a = 0 (we shall use the expression 
v'(0) = wi — (a/2)wo when needed for clarity). Moreover note that the kernel 
K{t) is square integrable. I 


We explicitly note that if(wo,wi) £ iJQ(f7)xL^(f7) (respectively (wq, Wi) £ 
L^(f7) X iL“^(f7)) then Fi(t) belongs to L^(0, T; L^(f7)) (respectively to 
L^(0, T; iL“^(f7))) and it depends continuously on (wo,Wi) in the specified 
spaces. 
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Using formula (TT^ we get the following Volterra integral equation of the 
second kind for w(<): 



(18) 


where * denotes the convolution. 



(19) 


and u(t) is the solution of (H)) when Uq = Wq and Ui = Wi — (o/2)wo and 
F(t) =0 (recall a = 0, see Remark O; 



Now we observe that the function 



belongs to C'([0, T]; i?i(U)) n C^i[0,T]-, L^(n)) if (wo,wi) e H^{n) x L^{n) 
and to C'([0,T];L2(f]))nC'i([0,r];iJ-i(U)) if (wo,Wi) e L'^{n) x H-^{n) and 


this is precisely the regularity of u(t). The fact that the solution of a Volterra 
integral equation retains the regularity of the affine term proves items [I] and [5] 
in Theorem [S] Now we prove the direct inequality, i.e. the statement [3l In this 
proof, T > 0 is fixed and we use M to denote constants (possibly depending on 
T) which are not the same at every occurrence. 

Remark 6 We shall use the following observation: the initial conditions enters 
into the affine term (and Wi is transformed to Wi — (a/2)wo) when we used 
MacCamy trick. If the original equation to be studied is 



V’(O) = I 

V^'(O) = T] 

■0 = 0 on . 


then 0 solves the Volterra integral equation dT51) with Fi = 0 and u{t) = 
R+{t)^ + A~^R-{t)r]. I 

In the proof of the statement [3] of Theorem[3]we assume f = 0, Wq G Rg (U), 
wi € L^(Q) and F G L^(0,T; L^(Q)). Let 



u(0) = uo = Wo G Ro(U) , u'(0) = ui = wi - ^wo G L^(U). 
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Picard iteration applied to (fTSl) gives 


w{t) 


u(t)+^ ^ f L{t — s)u{s) ds + A ^ 
Jo 


+ 0 O 

_^-n+2 

_n—2 



(t) (20) 


where denotes iterated convolution. 

We prove the direct inequality in the case wq G wi G 'D{Q) 

and F G V{n x (0,T)). In this case, w G (^([O,T];i?Q(fl) n H'^{Vl)) = 
(^([O, T]; dom A). The inequality is then extended to wq G Hq {Q), Wi G 
F G L^(0,T;L^(f2)) by continuity and density. 

Using dm) we see that there exists M = M(T) such that 


l’^’''^li2((0,T);L2(r)) <^(^l'7~u|^2((o_T);L2(r)) + 

+ T J L{t — s)u{s) ds^ 

+ D* 


+ 


\n—2 


L2(0,T;L2(r)) 
2 


L2(0,T;L2(r)) 


The known properties of the solution u of the elastic system © implies that 
(for a possibly different constant M = Mt) 


ITuI 


L2((0.T);L2(r)) 


+ 


'' +00 




< 


L2(0.T;L2(r)) 

< M (^|wo|^i(Q) + |wi||2(n) + |F|ii(o,T;L2(n)) 
In order to complete the proof we show that we have also 

rt 


T[A 


1-1 


J L{t — s)u(s) ds 


< 


L2(o,T;L2(r)) 


< M (|wo|^i(Q) + |wi|i2(n) + |F|ii(o,r;L2(Q))^ . 

f L{t—s)u{s) ds = A~^ [ L(t—s)u(s) ds+A”^ f L{t—s) f i?_(s—r)Fi(r) dr ds 
Jo Jo Jo Jo 


We note that 
rt 

A-^ 

lo Jo Jo Jo 

and the required inequality holds for the second addendum. The first addendum 
is 

pt 

A-^ 


= A-^ 
+ A -1 


/ L{t — s)u(s) ds = 

Jo 

' pt pt pt—r 

/ bR-{t — s)R^{s)uo ds / K{r) / R-{t — r — s)R^{s)uo dr ds 
Jo Jo Jo 

pt pt pt — r 

/ bR-{t — s)R-{s)ui ds / K{r) / i?_(i — r — s)i?_(s)ui dr ds 

Jo Jo Jo 
















Using (HI, we see that: 


■ pt pt—r 

TA~^ / — s)i?_(s)ui ds + / K{r) / R-{t — r — s)R-{s)ui dr ds 

Jo Jo Jo 

< M (|Wo|z,2(Q) + |wi|i 2 (Q)) < M ^|Wo|^fl(Q) + |wi|i 2 (Q)^ 

(recall m = wi — (a/2)wo). 

We consider the term containing Wq using 

R- {s)R+ (r) = i (i ?_{s + r) + R-{s - r)) . 


So we have 

TA~^ [ R-{t - s)R+{s)wo ds = ^TA~^R-(t)wo+ l:D*R-(t)wo . 

Jo 2 z 

The second addendum, as an element of C([0, T]; L^(r)), depends continuously 
on Wo S i?o (fl). 

The function is the solution of (H]) with f = 0, F = 0 and initial 

data uo = 0, ui = wq £ C so that 

\rA-^R-{t)wo\l2^o,T;LHr)) < ^|woli2(o) < . 

The remaining iterated integral is treated analogously. 

See [T7] for similar arguments and [T^] for a different proof (in the case of 
the wave equation). 


3 The proof of controllability 

In this section we prove the controllability result in Theorem As noted, the 
proof is in two steps. Let Re{T) = x L^(fl). Then we prove 

• Rv{T), at the same time T (and so also at larger times), is closed with 
finite codimension. 

• if e > 0 and </> T Rv(T+e) then </> = 0 and so Rv{T+e) = iL“^(n) xL^(n) 


3.1 The first step 

The first step is simple and it is based on the Picard formula (l20ll . Let wq = 0, 
wi = 0 and F = 0 so that u = u. We consider the operator f (-;■ w^(T) which, 
according to (EOl), is the sum of three terms, which are linear and continuous as 
functions off £ L^(0, T; L^(r)). The first one is 

f u(T) = u‘'(T) : 


this operator is surjective by assumptions (in fact, even more: f 
(T),u*^(T)^ is surjective). The third operator 


i^A-^ 


+ 00 

_n—2 



* U 


(T) 
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takes values in dom A = Hq{^1) n and so it is a compact operator from 

L 2 ( 0 ,T;L 2 (r)) to L^{n) . 

The operator D is compact and so the operators of f G L^(0, T; L^(r)) 


i 1 -^ A ^ f L{T — s)u^(s) ds = — f L{T — s) f i?_(s — r)Di{r) dr ds 
Jo Jo Jo 

with values in LJ{Vl) is compact. 

The function f h-)- w'{T), with values in is treated analogously. Note 

that 

—-A~^L(t)u = bR+(t)u+ / K(r)R^{t — r)u dr 

dt Jo 

belongs to C{[0,T]; L'^{n)) and L'^{Vt) is compactly embedded in and 

so the transformation 


— 

dt 


f L{t — s)u(s) ds 

Jo 


J t=T 


from L^(0, T; L^(r)) to H ^(fl) is compact. The same holds for the derivative 
of the last term in (EHl) computed for t = T, which is 





r) 


bg{r) + / K{t — r — s)g{s) ds 

Jo 


+ CXD 

g = U + ^ _4-n+2^(*(n-l)) ^ ^ ^ 
n—S 


dr , 


In conclusion, the reachable set Rv(T) is the image of the operator f i—>■ 
(T), u*^(T)^ , which is surjective in H~^{n)xL'^{n) since we assumed that 

system m is controllable at time T, perturbed by the addition of a compact 
operator. This implies that Rv{T) is closed with finite codimension. 


3.2 The second step 

The second step is more involved and requires several substeps. 

3.2.1 Step 2-1: the elements of 

We elaborate on formula (EOl) which can be written as 


+00 

w(t) = u(t) + / it(s)u(t — s) ds , it(t)v = ^ . (21) 

•^0 n=l 

We study the reachable set, hence we assume F = 0 and null initial condi¬ 
tions, so that 

u(t) = u(t) = .Av(t), v(t) = f R-{t — s)D {—f{s)) ds. 

Jo 

Of course, in the study of the reachable set, —f can be renamed f. 
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We characterize the elements (^, 77 ) G x L^(n) which annihilates 

Rv[T) (here T > 0 is arbitrary). I.e. we characterize those elements (4, 17 ) G 
Hq{Vl) X such that 

(((w*')'(T),|)) + f w^{x,T)r]{x) dx = 0 
Jn 

for every f G L^(0, T; L^(r)). The crochet denotes the pairing of and 

i/Q(r2). This set of annihilators is shortly denoted [Rv{T)]^ and we write 
itv)^Rv{T). 

Note that if it happens that w'(T) G L'^{Vl) then 


((w'('r),0)= f vf'{x,T)^{x) dx . 

Jn 

Of course we can study Rv{T)-^ by assuming / G ^(r x (0, T)) so that the 
following computations are justified. First we note that when / G ^(r x (0, T)) 
we have w(t) G L^(n) since 

rt rt rt-s 

w'(t) = ^ / ii_(s)Z3f'(t —s) ds + ^ / H{s) / i?_(r)Z3f'(t — s — r) dr ds . 

Jo Jo Jo 

Hence (in the last step we integrate by parts in time and we note that R-(t) 
and H[t) commute), 

((w'(r),0)= f ^{X)A r R.{T-s)Df{s)ds dx+ 

Jn Jo 


rT-s 


+ / ^(a;)^ / H{s) / R-{T — s — r)Df'(r) dr ds dx = 


f'{s) [D*AR-iT - s)|] ds dr+ 


r Jo 


dr dr = 


+ f f fir) D*A [ H{s)R-{T-r-s)^ds 
Jr Jo Jo 

= J Jo fix)i^D*A R+{T-r)^ + J^ H{s)R+{T - r - s)^ ds | dr dT . 

( 22 ) 

Note that the last integral makes sense thanks to the direct inequality, be¬ 
cause (l2T]l shows that the bracket is the solution of ([2]) with f = 0, F = 0 and 
initial condition w(0) = ^ G w'(0) = 0. 

Analogously, 

rT-s 


dx = 


[ ri{x) A f R-{T — r)Df{r) dr + A f H{s) f R-{r)Df{T — s — r)drds 
Jn Jo Jo Jo 

= f ( i{r)D*AR-{T — r)T] dr dT + f f i{r)D*A f H{T — r — s)R-{s)t] ds dT dr 
Jo Jr Jo Jr Jo 


J J f(r) a ^R-(T — r)T] + J H{T — r — s)A ^R-{s)'qds | 


dr dr. 


(23) 
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The bracket is the solution of @ with f = 0, F = 0 and initial condition 
w(0) = 0, w'(0) = r] G so that the brace belongs to L^(0,T; L^(r)) 

thanks to the direct inequality. 

We have _L i?y(T) when (1^ and (1^ sum to zero. Taking into 

account that the previous computation holds for every / £ T>iT x (0, T)) (which 
is dense in L^(0,T;L^(r)) we get 


D*A^^u{t)+ J — r)u(r) dr| = 0 , u{t) = + A ^R-{t)r]. 


Let 


■0(t) = u(t) + f H {t — r)u{r) dr . 

Jo 

The function 'ip{t) is the solution of the Volterra integral equation 

= u(t) + A~^ f L{t — s)^p{s) ds. 

Jo 


We compare with Remark [6] and we see that ^p{t) solves 


ij}" — Cij} + bij}{t)-\- f K(t — s)^p{s) ds , 
Jo 


^P{0)=^ 

1p'{0) = T] 

= 0 on dfl. 


(24) 


Note that the operator C is not in the memory term and that K{t) is square 
integrable. 

In conclusion, we can state: 

Theorem 7 We have T Rv{T) if and only if the solution '0(x,t) of the 

problem ^ 24 \ ) satisfy the condition 

D*Aip = 0. (25) 


As we noted, D*A = —T when applied to the elements of domA, in par¬ 
ticular when ^ G 'D{n) and rj G P(n). The direct inequality implies that 
T G C X L^(n),L^(0,T;L^(r)) and so T is the continuous extension of 

D*A to Hq{D,) X Lfip). Hence we can also state: 

Corollary 8 We have (|, rj) T Rv{T) if and only if the solution ip{x, t) of the 
problem j 24 l ) satisfies 

Ti/’ = 0 on T. 


3.2.2 Step 2-2: Fourier expansion and regularity of the elements of 

[Rv{T + e)]^ 

We recall that A is a selfadjoint operator with compact resolvent so that its 
spectrum is a sequence of eigenvalues and there exists an orthonormal basis 
{</>„} whose elements are eigenvectors of A. We denoted — the eigenvalue 
of (eigenvalues can be repeated, each one a finite number of times) and we 
assumed that the order has been chosen so to have A„ < A„+i. The bases 
of respectively iLo(fl) and iJ“^(f2) which correspond to the orthonormal basis 
{</>„} of the eigenvectors of A in L^(fl) are respectively {0„/A„} and {A„</>„}. 
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Let (^,77) G \RviT + e)]"'', (4, J7) G x L‘^{Vl). We can expand in series 

of the eigenvectors and we get 

+ 00 +00 

^ ^ <j>^{x)^ri 7 Tli^x) = ^ ^ 0^(x)?7n 5 {'^n^n} € ^ 7 {^n} ^ ^ ■ 

n=l n=l 

(26) 

First we prove the following result (which holds also with e = 0): 

Lemma 9 Assume Re{T) = x L^(f2) and e > 0. Let {^,ri) G 

[i?y(r + e)]-^ and {^,r]) 7 ^ 0. Then at least one of the series in Ii26]) is not 
a finite sum. 

Proof . Let us denote ipit) the solution of Eq. with ^p{0) = ^ and '0^(0) = rj 
and let us assume that both the expansions (I26p are finite sums, of N terms at 
most. We expand also the solution if{t), 

N 

0(0 = [0°(0?n + 0i^(O??n] 

n—1 

where both 0 °(t) and if\{t) solve the scalar equation 

0n(O = -^ni^nit) + hlfnit) + f K{t- s)0„(s) ds 

Jo 

with initial conditions respectively 

0 °( 0 ) = 1 , ( 0 °)'(O) = 0 , 0 i = 0 , ( 0 i)'(O) = 1 . 

The orthogonality condition is 

N 

0 = Tilj{t) = (Tcf^ix)) [0°(OC„ + 0 i)(O? 7 „] (27) 

n—1 

This sum cannot have only one nonzero addendum. To see this, let = 0 
and rjn — 0 for n ^ uq and either fno 0 0 or ijng 0 0 (or both). Computing 
either (EZl) or its derivative with t = 0 we get 

{rcf^,{x))=0. (28) 

It is a fact that the equality (1^51) cannot hold if L has been chosen in such a 
way that the elastic system m is controllable at some time T. The proof, based 
on the inverse inequality of Eq. o, is the same as that with d = 1 and can be 
found in [7] see also [HI Lemma 4.3]. 

Even more: the nonzero elements of j^Tp must correspond to at least two 
different eigenvalues. In fact, Eq. (EH) computed with t = 0 gives 



and if every 0„ correspond to the same eigenvalue —then ^n‘/’n(^) i® 

an eigenvector of A. As we noted, the equality is impossible. 
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Analogous argument if the derivative of both the sides of (l?7l) is computed 
with t = 0. 

Now we prove that if the sum is finite then we can reduce ourselves to have 
a sum of only one eigenvalue, and we proved that this is not possible. In fact, 
computing the second derivatives of (EH) we get 

N 

= 0 . (29) 

n —1 

In fact, the other terms appearing in the computation of the derivative sum to 
zero thanks to EZl). 

We subtract ((?7ll . multiplied with and (l29ll (as in [17]) so to obtain a 
new sum like (EZj) but with at most — 1 terms: 

A^-l 

{V'°(t) [A^ - A^] + V'i(t) [A^ - A^] r 7 „} = 0 

n —1 

After a finite number of iteration we remain with terms which correspond 
to the same eigenvalue (possibly, only one term) and we have seen that this is 
impossible. I 

Now we have this information, that at least one of the series EHj) is not a 
finite sum. So, the orthogonality condition (Ell) has to be replaced with 

+ 00 

0 = rV’(i) = + V'l^(i)» 7 n] (30) 

n—1 

(the exchange of T and the series is justified by the direct inequality). We are 
going to prove that also this case is impossible and so it must be ^ = 0 and 
T] = 0, but now we need the assumption that the purely elastic system o is 
controllable at time T and that (4, 17 ) € [RviT + e)]'*' with e > 0; i.e. that the 
equality (1501) holds in L‘^{0,T + e;L^(r)). The condition e > 0 is used in the 
proof of the following lemma: 

Lemma 10 Let Re(T) = x L‘^{Ll) and let e > 0. If {^,ri) € 

[RviT + e)]"*" then 

Cn = , Vn = ^ , {fn} S , {fjn} & l^ ■ (31) 

Proof . We expand in series of the eigenfunctions 4>^ the solution u(t) of the 
purely elastic problem El)-® when F = 0, Uq = 0 and Ui = 0. Using ® 
and ® (recall that —f was renamed f) we get: 

+ 00 

u(t) = ^ 0„(a;)w„(t), u'f = -X^^Un + 

n—1 


J^iTcfj-m dr. 


Hence, 


Unit) = ^ [ sinXnit-s) f (T</>„) ■ f(s) dr 
An Jo Ur 


d5 
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and we have the following expansions of u(t) and u'(f) (we compute with Ug = 0, 
Ui = 0 and F = 0): 


+ 00 


U(t) = -J2^n{x) 


n—1 
+ 00 


r</-„ 


sin A„(t — s) ) • f(s) dr ds , 


u'{t) =-^^{Xncj)„{x)) j J cos A„(t - s)^ • f(s) dr ds 

So, controllability is equivalent to the surjectivity of the map 


f 


A 

fTct>, 


Jo Jr \ 

• T 


COS A„s ) ■ f(r — s) dr ds 


n 


/O Jr \ 


sinA„s ) ■ f(T — s) dr ds 




Here l'^ = N = 1, 2, .... This transformation is continuous since f i—>■ 

(u'(T),u(T)) is continuous from L^(0, T; L^(r)) to x L^(f2). 

As usual with Fourier series, it is convenient to introduce 

A?.„ = -Xl , = 4>n 

and we see that controllability of the purely elastic system is equivalent to 
surjectivity of the following operator (here Z' = Z \ {0} and P = Z^(Z') = 


e £(L^ (0,T;L2(r)) . Mf = 


fTclj. 


V A, 


n S 


f(T - s) dr dsj 
(32) 

(see Lemmas 4.6 and 5.1 in |16]1: controllability of problem m and m is equiv¬ 
alent to the surjectivity of the bounded operator M in iS2[l from L^(0,T; L^(r)) 
to Z^(Z'; C). In turn, this is equivalent to the fact that the sequence 






where 

An 


is a Riesz sequence in L^(0,T; L^(r)) i.e. it can be transformed to an orthonor¬ 
mal sequence using a linear bounded and boundedly invertible transformation. 

We need the following pieces of information (see m ch. 3] for details on the 
Riesz sequences): 

Lemma 11 The following properties hold: 

1. if {cn} is a Riesz sequence in a Hilbert space H then '^anCn eonverges 
in H if and only if {an} G P; 

2. is a Riesz sequence in L'^{0,T; H) (H is a Hilbert space and 

kn & H) and if converges in L^(0,T -f e; H) to an H^ 

function then an = Sn/Xn and {i5„} € P. This result requires e > 0. 
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We continue the proof of Lemma[Tni we go back to examine the orthogonality 
condition =0 which can be written as 


+ 00 


^ = 0 , Zn{t) = Xu [tpnit)^n + tpnit)Vn 


( 33 ) 


n—1 


Hence Z„(f) solves 


Z'u — ~X^Z„ + bZu+ f K(t —s)Zu{s) ds , Z„(0) — , Z'^(0) — Xn^jn ■ 

Jo 

(34) 

So we have 

b /■* 

Z„{t) ={Xn^n) cos Xnt +r]n sin Xnt + — / sin A„(t - s)Z„(s) ds+ 

Jo 

+ — f sinA„(t —s) f K{s — r)Zn{r) dr ds (35) 

Xu Jo Jo 

(note that and ijn are real numbers and that {A„^„} G P because ^ G (H)). 

Gronwall inequality shows that the sequence of continuous functions {Z„(t)} 
is uniformly bounded on compact intervals. 

We have 


(Xn^n) cos Xui Tju sin Xui — e Cu “t e Cu ^ 
We introduce the notations 


Cn — 2 ■ 


Uu — c ^ Cu c ^ Cu j Su — sin Xut , Cu — COS Xui , Hu — bSu “t“ Su * H . 


iXfit 


Three steps of Picard iteration applied to (1551) give 

4 


Zu — Uu + T —Uu * 
An 


1 


V- vH. 


■-U) 


.l/=l 




where = M„(t) are continuous functions, and the sequences {M„(t)}, 
{{l/Xu)M^{t)} are bounded on bounded intervals. In fact, 


Mu = Uu* 


V 

M=5 


w ^3 =1= Zju ■ 


The orthogonality condition (|55)) takes the form 


+ 00 


+ 00 


- V = V '5'n^G, 

—^ ^ An 


^ ff*0 

i.= l 


+ 00 


J2^nj^Mu. ( 36 ) 


In fact, we can distribute the series since every one of the obtained series con¬ 
verges since 1) the last series converges because the estimate m implies that 


+ CXD 


1 


n=l " 
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2) the previous series converge since Sn *Un is a, linear combination of exponen¬ 
tials multiplied with a polynomial of degree at most 1 and our assumption 
on the time T implies that is a Riesz sequence. 

The right hand side of (1361) is an function on every interval. In fact, 
computing the derivative of the right hand side termwise we get the following 
sum of L^-convergent series, where Vn = U^ = (iA„c„) — (*A„c„): 


+ 00 ^ 

^ ^ ^nT (Cn ^n) 

1 ^7; 

n—1 
-\-oo 


V- vH. 

\7'-l ' 


*{v) 




Vn I * 


.17=1 

4 


V — 

117-1 


*(l/) 


+ 00 


/\,n 


(37) 


Convergence of the series is seen by an argument similar to the one used 
above, thanks to (HU and to the fact that {{1 /Xn)M!^{t)} is a bounded sequence. 
So, we can apply Lemma [TT] and we conclude 

Cn — T j G I hence Xn^n — T i — T ^ ^ ^ ^ ■ 

'^n 

(38) 

We equate the derivative of the left hand side of (l36ll and (1^ and we find 
the equality 


+ 00 


+ 00 




n—1 
+ 00 


/ 7 Ii7-l''''n 
17=1 


-b 




/ 7 117-1 
17=1 


+ 00 


/\,n 


(39) 


Using 

1 - 1 i ■ - 

Cn Cn — T (c-n, H“ C-n,) , — Ui (i) = "T [c ^ Cn 6 Ctj,] 

'^n 

we see differentiability of the first and second series on the right hand side. 
Using again the estimate m we see that also the last series is differentiable, 
since the explicit expression of (t) is 




(c„ + c„)y 

An 


fl—b 


+ ^ (C„ + * if) * 




-b 


In conclusion, the left hand side of dMl) is an function on every interval 
[0,T]. Applying again Lemma ITT] we get 


Cn - 1 Cyi , {Cyj} G I 

An 


SO 


that Xji^n — 2 (cn “b , 


_ “* U 
- 12 
71„ 


0 


as we wanted. I 


Note that we proved also this result: 

Theorem 12 The first and the second derivatives of the series in i33\) can be 
computed termwise. 
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3.2.3 Step 2-3: end of the proof 

We use Theorem [T2 the second derivative of (155)) computed termwise gives the 
equality 


+ 00 


+ 00 


= -J2^n - M{t- s) 


n—1 n—1 

so that we have also 
+ 00 


+ 00 


^ {XlZ^is)) 


ds = 0 


+ 00 




n—1 


n—1 


'^n 


= 0 . 


(40) 


Let Ni be the first index such that ^ 0. Combining (1551) and (1501) 

we get the new equality 


+ 00 


+ 00 


^ - A^) Z„(t) = ^ (A^^ - A^) Z„(t) 

n—Ni n—N 2 

+ 00 

n=N2 


X2 

A2 


- l) Vnipnit) 


= 0 


and N2 > Ni. So, the element (£ , 77 ^) S x L‘^{Vt) whose Fourier coeffi¬ 

cients are 


Cn = 


y A„ ’ 



'Hn , 


{AnCn} G , {77™} G 


is a second element in [i?y(T + e)]'*". 

The elements (^, rf) and are not colinear, since the Fourier coefficient 

of index Ni is nonzero for the first pair, while it is zero for the second one. 

Now we use the fact that (^i,J7i) S Hq{D,) x L^(n) and we iterate the 
procedure. We remove the first nonzero Fourier coefficient of rh) and so we 
get an element (^21 ^ 2 ) ^ ^d(^) ^ which belongs to [Rv{T + £)]"'■ and 

which is linearly independent from the previous ones. If the series of (^, rj) is not 
a finite sum then the procedure can be repeated and we get the contradictory 
statement that \Rv{T + e)]'*“ has infinite codimension. 

We combine this fact with Lemma |9] and we get ^ = 0, 77 = 0 as we wanted 
to prove. 
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